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Abstract: Single silicon vacancy (SiV) color centers in diamond 
have recently shown the ability for high brightness, narrow bandwidth, 
room temperature single photon emission. This work develops a model 
describing the three level population dynamics of single SiV centers in 
diamond nanocrystals on iridium surfaces including an intensity dependent 
de-shelving process. Furthermore, we investigate the brightness and pho- 
tostability of single centers and find maximum single photon rates of 6.2 
Mcps under continuous excitation. We investigate the collection efficiency 
of the fluorescence and estimate quantum efficiencies of the SiV centers. 
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1. Introduction 



Single color centers in diamond are auspicious for applications as solid state single photon 
sources (for a review see HI). Silicon vacancy (SiV) centers are especially promising due to 
their spectral properties such as the concentration of the fluorescence in a narrow zero-phonon- 
line (ZPL) with a room temperature width of down to 0.7 nm |2|. Furthermore, their emission 
is situated in the red spectral region at 740 nm, in a wavelength range where the background 
fluorescence of the diamond material is low 0. Studies using single SiV centers created via ion 
implantation in natural diamond, however, revealed a low brightness (approx. 1000 cps) lHISl. 
More recently, bright single SiV centers with up to 4.8 Mcps, created in situ, i.e., during the 
chemical vapor deposition (CVD) growth of randomly oriented nanodiamonds (NDs) and 
(001) oriented heteroepitaxial nanoislands (NIs) |6ll on iridium (Ir) films, have been observed. 
The origin of the enhanced brightness has not been fully explored up to now. In this context, it 
is crucial to investigate the collection efficiency obtained in the system as well as the quantum 
efficiency of the SiV centers. 

In this paper, we analyze several bright SiV centers in detail. The investigated SiV centers 
are hosted by NDs or NIs on Ir as introduced above (for sample details see ||2] |6|). Note that 
we additionally use a ND sample that has been grown with slightly modified CVD parameters 
compared to (55 min growth duration, 0.4% CH4) and contains slightly larger nanodia- 
monds (220 nm mean size). The SiV centers in this slightly modified sample showed similar 
characteristics. We extensively investigate the population dynamics of 14 single SiV centers 
including the saturation of the photoluminescence and explore the underlying level scheme, 
verifying a model including intensity dependent de-shelving which we suggested in 121 based 
only on the analysis of the population dynamics of one single SiV center Furthermore, we 
characterize the photostability of single SiV centers. We calculate the collection efficiency for 
the fluorescence of SiV centers on Ir. From this calculation and the maximum excited state pop- 
ulation, we estimate the quantum efficiency for the ZPL transition and indicate possible origins 
of non-radiative decay. 

2. Brightness of single SiV centers 

As an important figure of merit for a single photon source, we first determine the maximum 
photon count rate I^. For this experiment, we use a confocal microscope setup which is de- 
scribed in detail in l2]j6]|. For the randomly oriented NDs, excitation at 671 nm was employed; 
for the (001) NIs, an excitation wavelength of 695-696 nm was used (unless otherwise stated). 
As single SiV centers show preferential absorption of linearly polarized light f^, we employ 
the optimized linear polarization direction for excitation. All count rates are corrected for the 
dark counts of the setup. In the presence of background fluorescence, the fluorescence rate / of 
a single emitter as a function of the excitation power P is described by 

/ = /oo „ + ChackgrP- (1) 

" + "sat 

Using Eq. ([U, we fit the saturation curves for single SiV centers and obtain the saturation 
powers Psat and maximum photon rates /«, summarized in Fig. [T] We find a mean value for loc 
of 1 .5 ± 1 .4 Mcps in the randomly oriented NDs and of 1 .5 ± 2 Mcps in the (001) oriented NIs. 
The high standard deviation illustrates the variation in brightness of the emitters which will be 
discussed in detail in Sees. |3] and |5] The highest loo obtained from the fits is 6.2 Mcps. Thus, 
the single emitters observed here are the brightest color centers to date under continuous laser 
excitation. 

As apparent from Fig. [lib), also Psat displays a significant spread among different emitters. 
The highest value observed is 692 jiW, the lowest 14.3 /iW with a mean value of 105 ± 103 
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Fig. 1. Histograms of (a) maximum obtainable photon rates at saturation /oo (b) saturation 
powers Psai- Histograms take into account emitters from randomly oriented NDs and (001) 
NIs. NDs; excitation at 671 nm. NIs: excitation at 695-696 nm. 



/zW for the NDs and 175 ± 229 /iW for the NIs. By taking into account the transmission of 
the laser light through the microscope objective and the spot size of the focus (1/e^ radius) of 
approx. 0.5 /Jm, we estimate the intensity maximum of the focused beam impinging onto the 
color centers. The highest and lowest values of Pja, correspond to an intensity of 58.8 kW/cm^ 
and 1.2 kW/cm^. Thus, 2.1x10^-' photons-s^'-cm^^ impinge onto the emitter with the high- 
est Psat (excitation at 695 nm), whereas only 4.1 x 10^' photons s^' cm^^ are present for the 
emitter with the lowest Psm (excitation at 671 nm). 1 /iW at a wavelength of 671 nm (695 nm) 
coiTesponds to 2.87 (2.97)xl0^'' photons-s^'-cm^^. Note that the estimation of the intensity 
does not take into account the iridium surface, any losses (reflection, scattering, absorption) due 
to the nanodiamond or the occurrence of localized modes in the nanocrystals (discussion see 
below). The observed values of Psat are significantly smaller compared to previous studies using 
single SiV centers implanted into natural bulk diamonds: Wang et al. lH report P^at = 6.9 mW 
(excitation 685 nm, comparable or even tighter focussing of the excitation laser). The more ef- 
ficient excitation, first, might arise from local field enhancements at the site of the SiV centers: 
In (spherical) NDs with sizes comparable to the wavelength of the excitation/fluorescence light, 
resonant modes (Mie resonances) of the light field can develop [71. The excitation laser light is 
coupled into these modes; the resulting field distribution excites the color center Q. Depending 
on the position of the color center, it experiences a high or low excitation light intensity. The 
enhanced or reduced intensity, compared to the situation in bulk diamond where the Gaussian 
focus of the laser determines the local intensity, leads to a lower or higher Psat- This effect is 
indistinguishable from an altered absorption coefficient as it results in a change of P^at as well. 
Note that as the NDs/NIs used here are not spherical, the formalism used in Ref. [71 cannot be 
straightforwardly applied. Selecting emitters from a confocal fluorescence map introduces an 
experimental bias: one preferably selects bright single centers most probably experiencing high 
local intensity and thus efficient excitation. Thus, the histograms of Fig. [T] tend to summarize 
SiV centers which experience efficient excitation despite the fact that the sample might contain 
centers where local fields introduce less efficient excitation. 

A second cause for a reduced P^at might be an enhanced absorption coefficient. The ground 
state of the SiV center was reported about 2.05 eV below the conduction band edge |8l, ex- 
cluding an excitation of the color center's electrons into the conduction band for our 1.78 eV 
(695 nm) or 1.85 eV (671 nm) excitation. The excitation, therefore, most probably involves ex- 
cited vibrational states of the electronically excited state. Individual SiV centers show strongly 
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Fig. 2. g*-"' functions for varying excitation power. Consecutive g^^^ functions have been 
shifted for clarity (lowest excitation power g'^' function in each graph not shifted, higher 
(lower) excitation power g'^' functions shifted by 1 (0.5) each), (a) emitter ND3 (lower 
graph 0.03, 0.11, 0.15, 0.49P,fl,, upper graph 0.6, 0.8, 1.2, 1.7,3. 5,5. SPsa,). (b) emitter NDl, 
g(^) includes background connection (lower graph 0.08, 0.17, 0.28, 0.47/'5fl,, upper graph 
1.5, 3.3, 6.3, 11.9, 23.5, 32.7/',a/)-(c) emitter Nil (lower graph 0.0 1,0.02, 0.03, 0.07/'™,, up- 
per graph 0.2,0.3,0.6, 1.1, 1.4/',«,). 



varying vibronic sideband spectra in emission together with a varying overall emission into the 
sidebands Changes in the emission spectrum can also indicate changes in the absorption 
spectrum Q, thus potentially altering the absorption coefficient for a given excitation wave- 
length. In previous work by Wang m[5l, no sideband spectra were recorded precluding direct 
comparison. 

3. Intensity auto-correlation (g^^^) measurements 

Figure |2] exemplarily displays excitation power dependent functions for three individual 
SiV centers. The g^^^ functions have been normalized assuming that ^'^'(t) = 1 for long de- 
lay times T. All measurements reveal a distinct antibunching. Furthermore, the g^^^ functions 
exceed one for certain delay times (bunching). NI labels emitters located in nanoislands, ND 
labels emitters located in randomly oriented nanodiamonds throughout this work. For emitter 
Nil, a pronounced bunching already occurs at low excitation powers, while for emitter NDl it 
only becomes visible at elevated excitation powers. functions involving a bunching indi- 
cate a three level system. In a first approach, we use a simplified model depicted in Fig.|3]for 
the population dynamics: Levels 1 and 2 are coupled via a fast radiative transition (rate coef- 
ficient k2i), the photons emitted on this transition are detected to determine g'-^\ In contrast, 
level 3 acts as a shelving state populated via the rate coefficient k2i with the possibility of re- 
laxation into the ground state via ^31. As long as the emitter resides in state 3, no photons on 
the transition 2^-1 are detected. This simple model has been successfully applied to molecules 
involving shelving states BlOll . To obtain the g^^^ function, one solves the rate equations for the 
populations A^; resulting in 

= 1 - (1 +ae-l"l/"2 (2) 



intensity dependent 
de-shelving 




3 



Fig. 3. Schematic representation of the extended three level model employed to explain the 
population dynamics of single SiV centers, explanation see text. 



The parameters a, Tj and T2 are given by Q: 

Ti,2 = 2/(A±Va2-4B) (3) 

A = ^12 +^-21+^23 +^31 (4) 



B = ^12^23 +^12^31 +^21^31 +^23^31 (5) 
1 - 72^31 



^3i(T2-Ti) 



(6) 



The parameter Ti governs the antibunching, while T2 governs the bunching of the g^'^'' function. 
The parameter a determines how pronounced the bunching is. In contrast to Eq. (|2]i, the meas- 
ured g(2' functions display g'^' (0) / 0. For several emitters, this deviation is only due to the 
instrument response of the Hanbury Brown Twiss setup, i.e., in particular the timing jitter of 
the APDs (details see Q): Eq. (|2]i convoluted with the instrument response fully explains the 
measured data [see Fig.|2ta)+(c)]. For emitters ND3 and Nil, the deviation A^'2^(0) between 
the fitted value of ^'^'(0) and the measured datapoints is less than 0.05, witnessing very pure 
single photon emission with negligible background contribution. For other emitters, broadband 
background emission of the diamond material deteriorates the ^'^^ functions. For the spectral 
region of interest, the broad luminescence is attributed to sp^ bonded disordered carbon (in 
diamond films) introducing electronic states into the bandgap (e.g., ifTTl ) or to grain boundaries 
in the diamond material lfT2l . To take into account background luminescence, we follow Ref. 

ifTSi and include the probability pe that a detected photon stems from the single SiV center into 

(2) 

the fit of the measured correlation function (t) via 

gl;\z) = i + {g^^\z)-i)pl (7) 

Pe is obtained from the signal to background ratio in the saturation curves. From the fits of the 
g^^^ function, we obtain the excitation power dependent values of the parameters a, Ti and T2 
[see Eq. dU]. In the following, we aim at modeling the power dependence of these parameters 
and deduce the rate coefficients of the color center's level scheme. Examples of the measured 
power dependent parameters for six individual SiV centers are given in Fig.|4] 

In a first approach, we assume the rate coefficients k2i, ^23, ^31 to be constant, whereas ki2 
depends Unearly on the excitation power P: ki2 = aP. The assumption is justified, as the color 
center is excited to vibrationally excited states that typically relax within picoseconds to the 
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Fig. 4. Intensity dependence of thie parameters a, Ti and T2: (a) emitter NDl, (b) ND2, 
(c) ND3, (d) ND4 (exc. 705 nm), (e) Nil, (f) NI7. Blue solid lines represent fitting curves 
according to the intensity dependent de-shelving model. Red dashed lines represent the 
excitation power dependence from the model with constant rates. 



vibrational ground state in the excited state (state 2) lfT4l . Thus, the intermediate pumping levels 
do not accumulate population, therefore efficiently suppressing stimulated emission and thus 
saturation on the pumping transition. Under this assumption, A:2i, A:23, k-^ can be calculated from 
high and vanishing excitation power limiting values of a, X\ and T2 (details see ||2]), whereas o 
can be derived using the saturation power Psat 

{k2i+hl)Psa, ■ ^ 

As visible from Fig. ID the observed power dependence of a, T\ and T2 allows to estimate the 
high and low power limiting values of these parameters. Using the model with constant rate 
coefficients, we calculate the power dependent curves for a, Tj and T2, shown in Fig.|4]as red 
dashed lines. As evident from Fig. 21 the model with constant rate coefficients reasonably well 
describes the power dependence of Ti, and, for some emitters, the power dependence of a. 
Nevertheless, it totally fails explaining the power dependence of Xo'- '^2 significantly increases 
up to three orders of magnitude at low excitation power In contrast, the model predicts a nearly 
constant value of T2 at low excitation power. 



In Ref. 121, we tentatively suggested an extension of the simple three level model discussed 
above allowing to account for these deviations. We here verify this model using a larger number 
of emitters together with extended g'^' measurements. To extend the model, an additional in- 
tensity dependent transition process is included. Following approaches in the literature lfT5l[T6l . 
we assume that the process reactivating the color center from the shelving state (de-shelving 
process, rate coefficient k^i) is intensity dependent. In order to accurately describe our exper- 
imental results, however, we find that the simple linear excitation power dependence of ^31 
found in lfT5l[T6l has to be replaced by a saturation behavior; 

^31 = ;^ + 4, (9) 

where k^^ is an intensity independent part. As depicted in Fig. [3] the de-shelving might be 
realized via an excitation from the shelving state to higher lying states that returns the color 
center to the ground state (see lfT6l ). Such an excitation process might intrinsically exhibit a 
saturation behavior. 

For this new model, we calculate k23, k^i , k^^, d under the assumption ^21 + ^23 > ^31 : 

1^31 ~ ^ 

^23= 

^21 = i ~ ^23 

The superscript °° (") denotes the limit for high (vanishing) excitation power. To derive the 
power dependence of a, Ti and T2, two additional parameters have to be determined: c [see Eq. 
©J and (7. We can no longer obtain a from the saturation curve as it is no longer feasible to 
link the rate coefficients to f^r as in Eq. ([8]l. Instead, c and a have to be obtained from fits of the 
power dependent curves of a, Ti and T2: We employ Eqs. ([S])-®, together with the definitions 
of ki2 and A:3i [Eq. (|9]l], as fit functions with free parameters c and a. From the fits, we find that 
the power dependence of Ti is almost fully determined by a. Thus, we determine a from the fit 
of Ti (P). Using the resulting value of a as fixed parameter, we fit a{P) and obtain the value of c. 
To complete the description, we plot the resulting curve for T2(P) with these parameters. Using 
this procedure, we find the best accordance of the fits, displayed in Fig.|4]as solid blue lines, 
and the measured data. Table[T|summarizes the rate coefficients and fit parameters obtained for 
the individual SiV centers covered in Fig.|4] 
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Table 1. Rate coefficients deduced from the limiting values of a, Ti and T2 using the three 
level model including intensity dependent de-shelving and parameters c and ct obtained 
from the fits. For comparison also P,,,, is given. 
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As apparent from Fig. |4] using the extended model we obtain a much better concurrence of 
the fitted curves and the measured power dependence of a, Ti and T2: Especially for emitter NI7 



[Fig.gff)], ND2 [Fig.gtb)] and ND3 [Fig.gtc)] all curves are well described. For emitters ND4 
[Fig.Htd)] and NDl [Fig. HJa)], the rapid drop of T2 at intermediate powers is overestimated. 
This has been observed for several other emitters. Additionally, the power dependence of a can 
only be qualitatively described using the model with intensity dependent de-shelving for emitter 
ND4 [Fig.Ud)]. For Nil [Fig.Ue)], an extraordinary behavior of a is observed, including an 
increase of a at very low powers. Note that Nil is the brightest color center observed delivering 
/oo = 6.2 Mcps. 

In the following, we shortly summarize the sources of errors in the data evaluation and their 
consequences for the obtained rate coefficients kij. First, comprises a comparably large un- 
certainty as this rate coefficient is determined by T2 [see Eq. (fTOliI: The estimation of t| from 
the very steep curves at low excitation powers is challenging. Furthermore, a is often rather 
small at low excitation power (weak bunching), thus the proper determination of Tt is demand- 
ing. Second, k2i, in contrast, can be reliably assigned as it is determined by rj* [see Eq. (fTslll: 
tJ* can be precisely estimated as Ti (P) has a comparably small slope for low excitation powers. 
Third, A:23 includes a moderate uncertainty as it is governed by and a° [Eq. ([T2tl: Fig.Ua) 
illustrates the challenge in finding a°°. The power dependence of a includes two datasets: The 
one marked with filled squares (filled dots) has been obtained including (excluding) background 
correction for the g'^^ function fitting. Both fits well describe the measured g'^' function and 
yield very similar values for Ti and T2, respectively, but differing values for a: The instrument 
response washes out the g^^' function, thus a steep slope of g*^^-* in combination with a high 
absolute value close to zero delay leads to an increase of g^^-'(O) similar to the modification 
owing to background fluorescence. Thus, an uncertainty in the background correction induces 
an uncertainty in 0°°. In contrast, can be extracted reliably from the measured data as a clear 
convergence toward a constant value is observed for most emitters. 
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Fig. 5. Histograms of rate coefficients k2\, kri, ^31 and d obtained from the model of satu- 
rating de-shelving. The histograms take into account emitters from NIs as well as randomly 
oriented NDs. 



In summary, the measurements unambiguously reveal the presence of a shelving state 
and verify the existence of an intensity dependent de-shelving path. However, at present we 
cannot identify the nature of the shelving state and its energetic position nor the transition 
responsible for the de-shelving. The rate coefficients obtained from the intensity dependent de- 
shelving model for 14 emitters are summarized in Fig. [5] It is clear that, for all emitters, kn is 
significantly higher than the other rate coefficients. The significant spread of k2i might be due 
to the local environment as well as a varying quantum efficiency of the transition (for a detailed 
discussion see Sec. |5]). k^^ is lower than 1 MHz for the majority of the emitters, indicating a 



long lived shelving state, c/, representing the high power limit of the de-shelving rate coeffi- 
cient, is at least a factor of 4.6, for most emitters even an order of magnitude, larger than ^3 j . 
Comparing the parameter c, indicating the saturation power for the de-shelving process, with 
the measured saturation power Psat, we find that for several emitters c and Psat have a similar 
value. The values are summarized in Tab. [T] One thus might suspect, that the saturation of the 
de-shelving transition determines the saturation of the fluorescence of the SiV center The pa- 
rameter (7 gives the absorption cross section of the SiV centers (for details on the conversion 
from excitation power to intensity see Sec. |2]i. The values of a given in Tab. [T] correspond to 
absorption cross sections of 1.4^.2x 10^'^ cm^. Note that the excitation has been performed 
with optimized linear polarization to address the single transition dipole moment of the SiV 
center The absorption cross section for the nitrogen vacancy (NV) center under 532 nm excita- 
tion, averaged over the possible orientations of the transition dipole moments, has been recently 
determined to be approx. 1x10^'^ cm^ iflTl . The absorption cross section for the NE8 center, a 
nickel-nitrogen complex, has been determined as 1.7xl0^'^cm^ in |18 | for 687 nm excitation. 
The values for o determined for the SiV center here thus exceed the absorption cross section of 
the NV center by two orders of magnitude. To further clarify this issue, further investigations, 
e.g., using pulsed laser excitation as in ifTTl are desirable to determine the absorption cross sec- 
tion independently. For A:23 , a large spread is observed ranging from 137 MHz to 1 .5 MHz, thus 
we conclude that the coupling to the shelving state strongly varies among different emitters. 
It is also apparent from Fig. |5]as well as Tab. [T]that ^23 is comparable to or even larger than 
A:3i, even at high excitation powers. Thus, population accumulates in the shelving state. The 
influence of the shelving state will be further addressed in Sec.|5] 

4. Photostability of single SiV centers 

For single photon generation using optical excitation, the photostability of the emitter is crucial: 
Permanent or temporal loss of single photon emission, i.e., photobleaching or blinking, under 
optical excitation limits the applicability of the single photon source. For single color centers 
in diamond, photobleaching has been reported in the literature for single centers emitting in 
the near-infrared spectral region lfT9l . for single NV centers in NDs ll20l and for a center emit- 
ting at 736.8 nm |5 |. However, none of the publications discusses the origin of the permanent 
bleaching. In addition, fluorescence intermittence (blinking) is possible. For single color cen- 
ters in NDs, blinking has been observed in Ref. Il20l and is attributed to trapping and release 
of charges on the surfaces of the NDs. To analyze the photostability of single SiV centers, we 
obtain time traces of the fluorescence rate (see Fig.|6]l. Based on the observed fluorescence sta- 
bility, we can very roughly arrange the observed emitters into three classes as discussed below, 
whereas emitters may belong to class 2 and 3 simultaneously. 

Class 1: emitters with fully photostable emission 

An example for the fluorescence time trace of a photostable emitter (ND4) is given in Fig.|6ta). 
These emitters are photostable for excitation powers far above saturation: E.g., emitter NDl 
has been shown to be stable under excitation powers up to 32Psat- The typical observation time 
for the intensity dependent ^'^^ measurements discussed above is about one hour, thus these 
emitters have shown to be photostable for at least one hour under continuous laser excitation. 
Roughly 20-30% of the emitters investigated in detail show full photostability. As visible, e.g., 
in the lowermost graph of Fig. |6jb), the count rate of the stable emitters nevertheless exhibits 
slow variations (timescale > 5 min) due to spatial drifts of the emitter out of the laser focus 
which can be undone repositioning the emitter 
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Fig. 6. Fluorescence timetraces of single SiV centers (a) SiV center with fully photostable 
emission (emitter ND4, excited at 695 nm). (b) Emitter with destabilization at higher ex- 
citation power and permanent photobleaching (emitter Nil), (c) Emitter with longer time 
intervals of fluorescence intermittence (emitter NI6). The count rate of each emitter has 
been calculated in time windows of 100 ms for the lowest excitation power (50 ms for 
higher excitation powers). 



Class 2: emitters exhibiting fluorescence intermittence 

Figures |6lb) and (c) give time traces for emitters exhibiting only partially stable emission due 
to fluorescence intermittence. We find dark times ranging from several 100 ms up to 2 min. 
Figure |6|b) and (c) indicate a general trend: The probability for blinking events is apparently 
higher at elevated excitation power, thus one might suspect that the transition to the dark state 
is induced by the pump light. For excitation powers below or close to saturation, also for these 
partially stable emitters almost constant single photon emission can be obtained as the blinking 
events are rare. We point out that in Ref. Q an individual blinking color center with a ZPL 
wavelength compatible with SiV emission wavelengths ||6l is reported, evidencing the blinking 
of SiV centers in single crystal bulk diamond. 

Class 3: emitters for which permanent photobleaching occurred 

Figure |6jb) shows an emitter for which permanent photobleaching occurred at elevated power. 
'Permanent' here means, in general, that for waiting times of at least 10 min (partially without 
laser illumination) no recovery of the fluorescence has been detected. For emitter Nil [Fig. 
|6|b)], prior to the permanent bleaching event, blinking was observed with a trend to enhanced 
blinking activity for higher laser powers. However, we also found emitters that were bleached 
without any prior sign of fluorescence instability /intermittence preferably at higher excitation 
powers and after longer observations times (e.g., emitter N13 bleached at 13Pvo,, observation 
time 1 hour). We point out that in Ref. Q, permanent bleaching of the investigated blinking 
center after one week of observation is reported. 

4.1. Discussion of the observations 

Blinking of color centers can be due to photoionization of the color centers as the charge state 
after ionization can be non-radiative or emits at a wavelength beyond the preselected spectral 



window (here 730-750 nm). We here are not able to verify this option as the detection efficiency 
of our setup for the weak luminescence at 946 nm found to originate from an alternative charge 
state of the SiV center II2TI is too low. For NV centers in NDs with 5 nm size, blinking has been 
interpreted in terms of the capture of electrons in surface traps 1201 . However, simultaneously 
it was found that NV^ to NV" conversion is not responsible for blinking or bleaching as no 
NV" luminescence was observed ||2Q| . The authors of Ref. 1201 use the analogy of the optical 
excitation to an exciton formation to explain this behavior; As long as the electron of the exciton 
is captured, no fluorescence occurs. Other authors, in contrast, suspect that the lack of excess 
electrons needed to charge/recharge NV centers in small NDs leads to photobleaching after 
photoionization l22l . The observation of photostable SiV centers shows that the SiV complex 
is, in principle, photostable under red laser excitation, thus, modifications of the color center's 
local environment have to induce blinking or bleaching. For the SiV centers, the corresponding 
mechanisms are not clear However, due to the enhanced probability for the centers to undergo 
blinking transitions at elevated excitation powers, we suggest that the transition to the dark state 
is photoinduced as also found for NV centers in l20|. It might also be possible that the color 
center undergoes spontaneous transitions from its excited state to the dark state. Thus, with a 
higher excited state population the rate for a transition to the dark state is enhanced. 

We here choose the photon energy of the excitation laser sufficiently low so that the (spatially 
localized) electrons bound to the SiV center cannot be excited from the ground state of the SiV 
center to (spatially delocalized) conduction band states. Delocalized states may promote the 
capture of electrons by traps in the vicinity of the color centers: An electron in such delocalized 
states has a finite probability to be found at the spatial site of a trap. In addition to the primary 
excitation process, one might think of further excitation of the bound electrons from the excited 
state of the color center (excited state absorption). Also a simultaneous absorption of two pho- 
tons is a possible route to photoionization (c.f. two photon absorption processes for NV centers 
see, e.g., l23l ). In both cases, trapping of electrons promoted to the conduction band could 
subsequently induce fluorescence intermittence. The nature of the trapping states is not clear 
for the SiV centers investigated here. One might think of surface states as, e.g., in Ref. l20ll : 
however, as we use CVD diamonds, the surfaces of the diamonds should contain significantly 
less graphite and disordered carbon as compared to the surfaces of the detonation NDs used 
in I20I . A second possibility might be trapping of the electrons at other impurity atoms, e.g., 
substitutional nitrogen atoms 1241 . 

Additionally, it is not clear whether the mechanisms responsible for blinking and permanent 
photobleaching are identical. As very long blinking times have been observed, it is possible that 
the blinking mechanism is also responsible for the 'permanent' bleaching and that a recovery 
of the fluorescence after long waiting times is possible but has not been observed. This might 
especially happen if the laser is not able to free the electrons from their trapping states and if 
the spontaneous, possibly thermal, escape from these traps is very unlikely. The latter argument 
suggests trapping states deep within the band gap. The observation of photostable SiV centers is 
very promising for the application of single SiV centers as single photon sources. Using surface 
treatments as well as further control of the impurity content might help to enhance the fraction 
of fully photostable SiV centers. 

5. Quantum efficiency of single SiV centers 

This section deals with the maximum single photon rates loo observed for single SiV centers 
and the effects limiting this rate. In particular, we aim at deducing the quantum efficiency r\qe 



of the SiV centers. loo for continuous laser excitation is given by: 



loo = ride, riqek2lN2{P ^ °o) = T],/,, T],,, (14) 

^ + :^ 

N^{P — !• °°) is the maximum steady state population of the excited state. k2i, k23, k^^^ and d 
are the rate coefficients obtained from the intensity dependent de-shelving model. T]^/;,, is the 
detection efficiency of the experimental setup. It is the product of the collection efficiency rjcoih 
i.e., the probability to collect an emitted fluorescence photon, and the internal efficiency of the 
detection setup Tj'J"^, i.e., the probability to detect a collected photon. Taking into account the 
transmission/reflection of all optical components, as well as the APD detection efficiency and 
the efficiency of the multi-mode fiber coupling in the employed confocal microscope setup 
(details see |l2]|6l), we estimate the internal detection efficiency rj'J'^^ of our setup as 25%. rjqe 
is the quantum efficiency of the SiV center, i.e., the probability for a photon emission upon 
a transition from state 2 to 1 (see Fig. [3]). First, we determine the influence of the shelving 
state on loo. For an off-resonantly pumped two level system, assuming a very fast relaxation 
to state 2 after excitation, full population inversion Af~(P °°) = I can be obtained. For the 
emitters discussed here, we obtain A^~(f — > °a) as summarized in Tab.|2] As apparent from Tab. 
12] the influence of the shelving state on /oo differs for individual emitters: For emitter ND3, 
A'^(P oo) is only lowered by a factor of two compared to the two level case. On the other 
hand, for emitter NDl, A^^(P °°) is smaller by nearly an order of magnitude compared to 
the off-resonantly pumped two level system. As apparent from Tab. |2] the shelving rate ^23 is 
always much smaller compared to kji - However, due to a slow depopulation of the shelving 
state, for several SiV centers, the shelving state accumulates most of the population, leading to 
a significant loss of brightness compared to a two level system. 

Table 2. Rate coefficients kjj, maximum excited state population N^{P — s- 00), maximum 
photon rate /oo and quantum efficiency rfq^ for individual SiV centers. To calculate the quan- 
tum efficiency, we use a collection efficiency of 78% (28%) for a parallel (perpendicular) 
dipole, corresponding to an emitter distance of 75 nm. 
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5.7. Calculation of collection efficiency and quantum efficiency: influence of the Ir substrate 

The radiation properties of SiV color centers in NDs/NIs on Ir are investigated by considering 
a point-like oscillating dipole near a metal surface 1251 . The orientation of the individual color 
center dipoles is unknown, we thus investigate the limiting cases of a dipole perpendicular and 
parallel to the interface in our simulations. Determining the exact position of the emitting SiV 
centers and thus their distance from the Ir surface is, in principle, not possible as they are cre- 
ated at an unknown instant of time during the CVD growth of the nanodiamonds, i.e., at an 
unknown position inside the CVD nanocrystal. For an estimation of the characteristic distance 
of an emitter from the surface, we assume the emitter to be located roughly in the center of the 
NDs/NIs, i.e., approx. 40-100 nm above the metal film. These assumptions seem to be crude 
approximations as we are neglecting the fact that the color center is inside a dielectric nanopar- 
ticle, which may affect the photophysics of the SiV center through changes in the spontaneous 
emission rate 1261 and different regimes of interaction with the metal surface 1271 . Nonetheless, 
since we are interested in a qualitative understanding of the role of the Ir surface and we are 
concerned with distances above roughly 50 nm, it turns out that our approach is sufficiently 
sophisticated to describe the most important involved phenomena. 

The dipole emits radiation in vacuo at A = 740 nm, where the dielectric function of Ir takes 
the value eji = — 1 8 + 25i ||28l . In practice, the Ir surface modifies the quantum yield by chang- 
ing the spontaneous emission rate and absorbing a fraction of the emitted light. If yo and Tjo are 
the intrinsic radiative decay rate and quantum yield of the SiV center respectively, an expression 
for the effective quantum yield reads 



(l-f?o)ro/7r + W'?a' 

Yr represents the modified radiative decay rate and T7a accounts for the rate /m- of energy dissipa- 
tion in the metal, Tja = Yr/ (7r + Ym)- The effective quantum yield for a parallel and a perpendic- 
ular dipole is shown in Fig.|7ta) for tjq = 5% (according to ||29l ) as a function of distance from 
the Ir surface. Note that the competition between Tja and Yv/Yo may lead to an effective quantum 
yield larger than tjo at certain distances if the dipole is parallel to the Ir surface. The effect is, 
however, modest being Yr at most a factor of two larger than 70. For very short distances, the 
increase of Ym for both dipole orientations due to near-field energy transfer, gives rise to the 
well-known phenomenon of fluorescence quenching |25i]. 

The collection efficiency is given by the fraction of power radiated in the solid angle deter- 
mined by the numerical aperture (NA) of the microscope objective, divided by the total radiated 
power which, in the presence of a metal surface, is limited to the upper half space. In short, the 
calculations are performed by expanding the dipole field in plane waves. Each partial wave 
fulfills the boundary conditions at the interface through Fresnel coefficients. Further details can 
be found in l30l[3Tl . The collection efficiency strongly depends on the radiation pattern, which 
is significantly modified by the Ir surface. An example for parallel and perpendicular dipoles 
located 80 nm above the metal is shown in Fig.|2tb), where the pattern for the respective dipoles 
in free space is added for comparison. In practice, the Ir surface channels the emission toward 
smaller angles, thus increasing the fraction of emission that can be collected by the microscope 
objective. By changing the distance, the radiation pattern varies differently for the two relevant 
dipole orientations, as it can be inferred from Fig.|2tc), where the collection efficiency using a 
microscope objective with NA=0.8 is plotted as a function of distance. It is found that a par- 
allel dipole is a favorable configuration for obtaining large count rates, where the fraction of 
collected photons can exceed 70% with standard optics for microscopy. We recall that for the 
case of SiV centers in bulk the large refractive index of diamond leads to poor collection effi- 
ciencies of only up to a few percent, for a dipole parallel to the interface (not shown). For the 
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Fig. 7. Emissions characteristics of a point-like oscillating dipole in vacuum near an Ir 
surface for orientation parallel (blue solid curves) and perpendicular (red dashed curves) to 
the interface. The emission wavelength is 740 nm and the intrinsic quantum yield is 5%. (a) 
Effective quantum yield as a function of distance from the Ir surface, (b) Radiation pattern 
for a dipole located 80 nm above the Ir surface. The thin curves refer to a dipole in free 
space, (c) Collection efficiency as a function of distance for a microscope objective with 
NA=0.8. 

typical emitter distances of 40-100 nm, we find a collection efficiency of approx. 75% (30%) 
for parallel (perpendicular) dipoles, each value varying within less than 10% [see Fig.|3c)]. 

In the following, we compare our experiments using NDs/Nls on Ir to previous experiments 
using SiV centers in bulk diamond in ||4]|5|. For a parallel dipole, the collection efficiency 
in the NDs/NIs on Ir is enhanced by approx. a factor of 20 compared to SiV centers in bulk 
diamond. On the other hand, the maximum photon rates !„ we find are about three orders of 
magnitude higher than observed in ||4]|5]. 7]f^j in Ref. JS] is 19% and thus close to the value 
of our setup (25%). However, it should be noted that fsl does not include a discussion of the 
transition rate coefficients similar to the discussion given here. The quantum efficiency of a 
single SiV center modeled as a two level system is estimated to be 0.5% in Q. Therefore, we 
tentatively suggest that the high brightness of some of the SiV centers in CVD NDs/NIs cannot 
fully be attributed to an enhanced collection efficiency but is also correlated to a slightly higher 
quantum efficiency for the in situ produced SiV centers as discussed below. 

5.2. Experimental estimation of the quantum efficiency 

We now use rjcoii calculated for a dipole located 75 nm above the Ir surface to estimate rj^g of 

the SiV centers according to Eq. ( fT4l i. We find values ranging from rjlje = 0.3% to 7]^ = 9.2% 
assuming the two limiting cases of parallel and perpendicular dipole orientations (Tab.|2]). Thus, 
the observed quantum efficiency is comparable to previous measurements on SiV ensembles in 
polycrystalline films yielding rji/e = 5% 1291 . The values determined here, however, do not 
straightforwardly represent the internal quantum efficiency, i.e., the probability for a radiative 
decay of the SiV center in bulk diamond. First, for the saturation measurements, from which we 
obtain I„c, only the fluorescence in a spectral window 730-750 nm is recorded, thus neglecting 
roughly 30% of the fluorescence emitted into sidebands and additional electronic transitions in 
the near infrared spectral range (see ll32l ). Thus, these transitions are erroneously considered as 
non-radiative and the quantum efficiency rjge will be underestimated by approx. 30%. However, 
as the fractional intensity of the red-shifted emission significantly varies for individual emitters 
||6| also the error for the estimation of rji^g varies. In addition to influencing ?7co//, the presence of 
the metal can also quench the fluorescence for emitters close to the surface, thus reducing rjgg. 
However, to estimate this influence, the unknown intrinsic quantum yield of the SiV centers as 
well as the exact distance to the metal has to be taken into account. Thus, t]^^ estimated here is 



the probability for a photon emission in a restricted spectral range for an SiV center above the 
metal surface. Furthermore, in addition to our simple model of a dipole in air above the metal 
surface, the ND can strongly modify the radiation pattern of the emitting dipole: As discussed 
in Refs. ||2l|33] and mentioned before, in spherical particles with sizes comparable to the wave- 
length of the fluorescence and excitation light, resonant modes (Mie-resonances) can develop, 
strongly modifying the radiation pattern and thus rjcoii- Ref- Ell demonstrates a variation of 
77^0;; between approx. 1% and 20% (spherical NDs on a sapphire substrate, size varies from 
50-200 nm). However, spherical NDs still simplify the problem as the nanocrystals mostly re- 
semble cubo-octahedral shapes and their exact size is unknown ||2]|6l. We cannot estimate the 
modification of Tj^,// due to possible resonant modes in our NDs/NIs of unknown shape and 
size. Therefore, we cannot state whether this effect leads to an overestimation or underestima- 
tion of rjge- In addition to these considerations, we emphasize that the rate coefficients and thus 
the excited state population include uncertainties (discussed above). Furthermore, the measure- 
ment of the maximum obtainable photon rates includes an uncertainty due to the contribution 
of background fluorescence. A quantum efficiency strongly deviating from 100%, in principle, 
might have several reasons: First, quenching of color center luminescence by the proximity to 
defect rich crystal areas has been reported in the literature: Refs. Il34l . If33]| and 1351 indicate 
quenching due to graphite and disordered carbon on ND surfaces, crystal damage as a conse- 
quence of heavy ion irradiation and structural defects or non-diamond carbon phases. The actual 
process leading to the quenching is not discussed or identified. Despite a high crystal quality 
of the NDs/NIs employed in this work, defects like dislocations are present and might possibly 
induce a quenching. Moreover, in a solid state host radiative transitions can be quenched by 
the direct emission of phonons (multi-phonon relaxation). The 1.68 eV transition of the SiV 
center equals 10.2 phonon energies (165 meV in diamond). Ref ||36l summarizes evidence for 
multi-phonon quenching for luminescent transitions in diamond with similar energy. An addi- 
tional quenching mechanism for color centers has been introduced by Dexter, Klick and Russel 
in 1371 . The non-radiative process is induced at a crossing point in the configuration coordinate 
diagram where the energy of a low lying vibrational state in the electronically excited state 
matches the energy of a highly excited vibrational state in the ground state and the color center 
can relax to the ground state without emission of radiation. Such processes might exist for SiV 
centers in diamond but have not been considered so far. 

6. Conclusion 

Single SiV centers in NDs and NIs on Ir films have been shown to exhibit high brightness 
under continuous excitation. We have developed a model accurately describing the three level 
population dynamics of these centers including an intensity dependent de-shelving process. 
SiV centers have been observed to retain photostability for excitation well above saturation, 
however, also blinking or bleaching centers have been found. The employed material system, 
nanocrystals on Ir, enables a high fluorescence collection efficiency exceeding 70%. With this 
observations, we estimate quantum efficiencies for single SiV centers of up to 9%. 
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